clear;close all;
dimX=1000;
repeats=1000;
%generate a matrix
x=rand(dimX,'double');

%calculate
tic;
testfun1(x,repeats);
toc;
